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^; Abstract 

^ ' We discuss effective field theory treatments of tlie problem of three particles interacting via 

^^ ■ short-range forces (range R <^ 02, with 02 the two-body scattering length). We show that forming 

a once-subtracted scattering equation yields a scattering amplitude whose low-momentum part is 

renormalization-group invariant up to corrections of ©(ii^/al)- Since corrections of 0{R/a2) and 

^ \ 0{R^ /a'2) can be straightforwardly included in the integral equation's kernel, a unique solution 

^ ' for 1+2 scattering phase shifts and three-body bound-state energies can be obtained up to this 

accuracy. We use our equation to calculate the correlation between the binding energies of Helium-4 

trimers and the atom-dimer scattering length. Our results are in excellent agreement with the recent 

three-dimensional Faddeev calculations of Roudnev and collaborators that used phenomenological 

inter-atomic potentials. 
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I. INTRODUCTION 

Zero-energy scattering between two particles can be characterized by one number: the 
scattering length 02- The zero-energy cross-section for scattering of two non-relativistic 
quantum-mechanical particles is: 

a = Anal- (1) 

This formula governs the cross section at center-of-mass momenta k which are significantly 
smaller than 1/R, where R is the range of the underlying two-particle potential. In the 
domain where kR -C 1 an effective field theory (EFT) with contact interactions alone can 
be used to organize the corrections to Eq. (jT)), and relate a to other observables, via a 
controlled expansion in a small parameter or small parameters. The nature of this expansion 
is, however, quite different depending on the relative sizes of the scattering length 02 and 
the range R. In the first case, 02 ~ R, one can calculate corrections to Eq. (jT)) via a 
perturbative expansion in 02 k. But, if 02 3> R — as is true in nuclear physics and a number 
of other quantum-mechanical systems — the expansion in ka2 breaks down very quickly. It 
is then more efficient to develop an expansion in the small parameters R/a2 and k R. Key 
to this modified expansion is that, if 02 > 0, the large scattering length is associated with a 
shallow bound state with binding energy 

B2 = ^. (2) 

(Here and in the rest of this work we set ^ = 1.) The existence of this bound state signals the 
onset of non-perturbative physics in the EFT. In order to describe physics at fc ~ 1/02 <C 
1/R we need an EFT expansion that already at its leading order generates the two-body 
bound state with energy B2 (or its analog on the second Riemann sheet if 02 < 0). The 
properties of this bound state can then be accessed by scattering particles from the bound 
state, i.e. examining the three-body scattering problem which we will refer to from now on 
as "1-1-2 scattering". It is this type of scattering which we will focus on in this paper. 

The most extreme example of the situation 02 ^ -R occurs when we take the limit 
02 -^ cxD, or, equivalently R —* 0. When zero-range interactions are involved the three-body 
problem is ill-defined without the introduction of an additional length scale, which can be 
chosen to be the scattering length for the 1-1-2 process, 03 [l|, |2, |3|. This means that in the 
limit a2 —>■ 00 the three-body system is characterized solely by the ratio a^/a2- Universal 
predictions for the low-energy (i.e. kR <^ 1) behavior of any system with a particular value 
of that ratio can then be generated — either using an EFT or quantum-mechanical models. 
Since the EFT generates predictions for these systems which do not depend on the details of 
the short-distance physics it provides a powerful way of unifying apparently diverse physics 
problems: dimers, trimers and tetramers of IIelium-4 atoms; deuterium, tritium, & Helium 
nuclei; and nuclei with neutron halos |j, y, la lD, |3| . These systems differ in size and spin 
structure but share the same essential feature of a large two-body scattering length. For a 
recent review of this physics, see Ref. |9|]. 

In reality, no system can exist exactly in the universality limit a2 -^ 00. Being able 
to calculate the deviations from, and approach to, universality is important for all such 
problems. That approach is controlled by the parameter R/a2- A first assessment of the 
impact that taking R/a2 7^ has on universality predictions was carried out by Efimov 10 , 
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ll|, and these calculations were recently systematized using EFT by Hammer & Mehen 



and Bedaque et al. |^, both of whom showed that the 0{R/a2) corrections to universality 



came from the physics of two-body scattering alone. This means that, up to corrections 
of 0{R^/a'^, the low-energy properties of the three-body system are determined by two 
numbers: the ratios R/a2 and 03/02- 

In this paper we take this argument one step further, and use EFT to show that the prop- 
erties of the three-body system are actually determined by those two ratios up to corrections 
ofO{R^/a\). In order to do this we develop the integral equation describing three-body scat- 
tering up to an arbitrary order in R/a2. It has been known for forty years that to zeroth 
order in this parameter, i.e. when 02 — > cxd, this integral equation is ill-defined [2]. In Sec- 
tion |n] we discuss the nature of this problem in the language of the renormalization group, 
and reiterate the argument of Ref. jj in order to show that the equation can be made 
well-defined by one subtraction, which depends on the three-body parameter 03. Therefore 
once 02 and 03 are known the integral equation may be used to generate the predictions 
for 1+2 scattering phase shifts, three-body bound-state energies, recombination rates, etc. 
mentioned above. Then, in Section UTTl we show that this same subtraction also makes the 
equation well-defined when 0{R/a2) and 0{R^ /a^ corrections are included in the kernel of 
the integral equation. This suffices to prove that only 02, 03, and R are needed to generate 
predictions for three-body system observables which are accurate up to 0{R? /a^). Bedaque 
et al. |J| and GrieBhammer |l4| have claimed that an additional three-body datum is neces- 
sary to achieve this accuracy, but our results show that renormalization does not necessitate 
the inclusion of this information in the calculation at 0{R^ /a^. 

These results are true for all of the physical systems listed above — indeed, for any system 
for which R/a2 is a small parameter. In order to demonstrate their utility we focus on 
an apparently straightforward application: to dimers and trimers of Helium-4 atoms. In 
Section IIVI we use our formalism to generate predictions for the binding energies of trimers 
of Helium-4 atoms. Since the atom-atom scattering length is 02 = (l04l^g)A while the 
typical van der Waal's forces between the Helium atoms have a range i? Ri 7 A, in this case 
the expansion parameter R/oo ~ 10%. The '^He trimer, tetramer, and several larger ^He 
clusters have been observed |l5l.ll6l|. but unfortunately there are no measurements of trimer 
binding energies, and so we use precise few-body calculations |l7|, |l8|, |l9|, |20|, |21[ , based on 
quite complicated inter-atomic potentials J22, |23| , to provide the two input ratios for our 
EFT. We then compare the EFT's predictions with the corresponding results obtained in 
Refs. jl9|, |20|. In fact, the numbers obtained for the atom-dimer scattering length in these 
two papers differ by more than the stated systematic error of the computations. Our EFT 



calculation confirms the result of Ref. 1£|, but we are unable to reproduce the result for 
the scattering length obtained in Ref. 20|. Indeed, the discrepancy with their number is 
two orders of magnitude larger than the expected R^/a2 error of our calculation. From the 
perspective of the EFT this renders the results of Ref. J20| suspect, since they are not in 
accord with the model-independent pattern of R/a2 corrections that we have computed here. 
With the addition of spin and isospin degrees of freedom, the formalism developed here 
can be applied to neutron-deuteron and (with the inclusion of Coulomb corrections) proton- 
deuteron scattering at low energies. In fact, the power counting used here for atomic systems 
has been quite thoroughly tested in the nuclear context. Calculations to NNLO (and several 
orders beyond in the A^A^ system) converge as expected |^, |2J, |25|. We will pursue the use 
of our formalism for such systems, and also develop a perturbative expansion of the solution 
of our integral equation, in future work. 



II. THREE-BODY SCATTERING AT LEADING ORDER 

In this section we will give a short introduction to the application of the EFT with contact 
interactions alone to the three-body sector and the renormalization of the integral equation 
that results after one subtraction. 

At very low energies all short-range interactions appear point-like and one can describe 
the physics of a system that contains only such interactions using an EFT built up from 
contact terms alone. The most general Lagrangian describing a non-relativistic system of 
identical bosons with short-range interactions is given by 

£ = ^t L + ^"J ^ - ^ii^^i^r - ^ii^w + • • • , (3) 

where the ellipses denote interactions with more derivatives and/or more fields. For the 
application to the three-body sector it is useful to rewrite this Lagrangian in terms of a 
particle field and a dimeron field, T, which has the quantum numbers of a two-body bound 
state 

C = tp^ [idt + — ]ij + AT^T - -^{T^ij + h.c.) + hT^Tij^ .... (4) 

2m I y 2 



It can be demonstrated J2a, |22| that these terms give the leading-order (LO) t-matrix for 
two-to-two scattering, in the context of the EFT expansion in powers of R/a2 developed in 
Refs. |28l . I29I |30|, |3l| . In the normalization we use in this paper that LO t-matrix is: 



Tim? E + B2 E + B2 



r^'\E) = —- J^J^^ ^—^, (5) 



where we have adjusted the constants g and A such that the two-body binding energy B2 
is correctly reproduced, and: 

S(^\E) = ^U + ^/^^^] . (6) 

Here, and in what follows, the superscript ^^^ indicates that this is a leading order result. 
Also 7 = ^JraB^ is the "binding momentum" of the two-body system. Equation (j2)) implies 
that at LO 7 = -. 

When embedded in the three-body system the two-body amplitude (0) generates the 
following equation for the 1-1-2 scattering amplitude, X^, provided that the three-body force 
h is, for the time being, ignored: 

Xf\q, q'- E) = 2Ze{q, q'; E) + ^ Aq"q"' 2Z,{q. q" , E) f'\^f,:'^"\ X,{q" . q'; E) , (7) 

Jo E-Y— + B2 

where £ denotes the relative angular momentum between the dimer and the third particle 

and5(0)(E;g") = ^W(i?-^). 

Considering only s-waves and defining 

07 ( ' T?\ -71 ' j^\ ^ ^ f q^ + q'"^ + qq' - mE \ 

2Zo{q,q;E) = Z{q,q;E) = ; log ^ 1 /2 'f ^ ' W 

qq \q^ + q'^ — qq' — mE J 



we reformulate the integral equation in Eq. ((Tj) via a principal-value integral 

K^'\q,q'-E) = Z{q,q'-E)+V [\q" q"'Z{q,q" , E)-^^^^^^^K('\q" ,q'- E) . (9) 

Jo E + B2- jq"^ 

Below breakup it is convenient to use this f^-matrix approach. (From now on we drop the 
subscript i, in order to avoid notational clutter. For the rest of the paper we will focus only 
on £ = 0.) The ii'-matrix is normalized such that it gives on-shell {E = ^ — B): 

ir(o)(fc,A;;E) = -^tan5. (10) 

87A; 

Using the boson-dimer effective-range expansion for the s-wave phase shift 

kcot5 = + -A;V3 + ... , (11) 

ag 2 

with as and r^ being respectively, the 1-1-2 scattering length and effective range, we find 

ir(°)(0,0;-i?,) = ^. (12) 

07 

The solution to the integral equation is strongly cutoff dependent. This is a manifestation 
of the fact that Eq. ((Tj) is not properly renormalized. Bedaque, Hammer, and van Kolck [3| 
have shown that one way to remedy this problem is to insert the dimer-boson coupling h 
of Eq. (j3)) into the kernel of the leading-order equation. This makes it possible to fix the 
behavior of the amplitude for arbitrary cutoffs by adjusting h such that a chosen three-body 
observable is reproduced correctly. 

An equivalent procedure is to use subtracted equations. In this approach we trade the 
three-body force h for a three-body observable which appears explicitly in the equation. 
The result is a formulation of the three-body problem with short-range forces in which only 
renormalized quantities appear. 

We begin by considering: 

4m /"^ 
K^'\q, 0; -B,) = Z{q, 0; -B,) -—J dq" Z{q, q"; -B,)S^'\-B,, ; q")K^'\q", 0; -B^) , 

(13) 
We can insert the required three-body input by demanding that the scattering amplitude 
reproduces the correct particle-dimer scattering length at threshold. The on-shell amplitude 
a.t E = —B2 is given by Eq. (fT^ and it obeys 



Am r 
K(^){0,0;-B,) = ZiO,0;-B,)-^ / dq" ZiO,q";-B,)S^'\-B„;q")K^'\q",0;~B,) , 

•J Jo 

(14) 

where 03 denotes the atom-dimer scattering length. By subtracting Eq. p4|) from Eq. (fT^ 



and using Eq. (fT^ we obtain the subtracted equation at threshold^: 
K^'\q,0;-B2) = ^-^ + A[Z]iq,0;-B,) 

Am f^ 

-^ / dq"A[Z]{q,q"--B,)S'^'\-B,;q")K^'\q",0;-B,) , (15) 
<J Jo 

where 

A[Z]{q, q'; E) = Z{q, q'; E) - Z(0, g'; E) . (16) 

This equation was first derived in a different notation by Hammer and Mehen J32]. As they 
showed, it is renormahzation-group invariant up to corrections which scale ~ q^ — a result 
that is in accord with the more recent analysis of Ref . J33| . 

Now we can determine the full-off-shell amplitude at threshold by exploiting the symme- 
tries of K. Before subtraction the full off-shell amplitude at threshold satisfies 

Am r^ 
K^'\q, q'- -B,) = Z(g, q'; -B,) -^ j^ dg" Z(g, q"- -B,)S^'\-B,- q")K^'\q\ g'; -B,) 

(17) 
Since Z{q,q';E) = Z{q',q;E), we also have K^^\0,q; -B2) = K^°'^{q,0; -B2). Using this 
fact as input, together with the solution of Eq. (fT3]) . allows the subtracted version of Eq. (|T7jl : 

K^'\q,q'--B2) = K'^'^\0,q'; -B^) + A[Z]{q,q'- -B^) 

Am r^ 
-— j dq"A[Z]{q,q"--B,)S^'\-B,-q")K^'\q\^--B,) ,(18) 

to determine the fully-off-shell amplitude at threshold. 

With Eq. (jTpi in hand resolvent identities may be used to obtain the subtracted amplitude 
at any energy |l3lj |: 

ir(0)(g,A;;E)=ir(0)(g,fc;-fi2) + 5W(g,fc;-52)+ / dq'q'^Y^''\q,q'-E)K^''\q',k-E) , 

Jo 

(19) 

where the second inhomogeneous term is given by 

B^'\q,k-E) = 5[Z]{q,k-E)+ [ dq"q"^K^'\q,q"--B,)r^'\~B2-.q")5[Z]{q\k-E), (20) 

Jo 

with 

5[Z]{q, q'- E) = Z{q, q'; E) - Z{q, q'; -B^) , (21) 

and t{E] q) = t{E - |^). The kernel in Eq.dTHD is 

r(°)(g, q'; E) = K^'\q, q'; -i?2)5[rW](E; q') + 6[Z]{q, q'; E)t^'\E- q') 

+ r dg" q"' K^'\q, q"; -B,)r^'^ {-B,- q")5[Z]{q'\ g'; E)r^'\E- q') 
Jo 

= K^'\q,q'--B2)6[T^%E-q') + B^'\q,q'-E)T^'\E-q') , (22) 



^ Note that we would still have obtained Eq. H15|l if we had started out with a scattering equation containing 
the momentum-independent three-body force h and considered the equations corresponding to (|13|l and 
lfn|) in that case. 



with: 

S[r^'mE;q')=T'~'\E;q')-r^'\0-q'). (23) 

The asymptotic behavior of the elements of Eqs. (fT^ -()23 p is as follows: 

r(°)(E;g')~^ for g' > fc (24) 

(5[r(°)](E;g')~-^ for g' > ^ (25) 

ir(o)(g,g';-52)~^ for g, g' » 7 (26) 

qq' 

K^'\q,q';-B,)r^— for g' ~ 7, g » g', 7 (27) 

g7 

5[Z]{q,q';E)r^-^ for g~g'»fc,7 (28) 

(In Eqs. ()24p - ()28p we omit factors of m since they are simply overall factors in the final 
answers obtained for these non-relativistic systems.) In the case of K^'^'' the behavior listed 
is modulated by an oscillatory function, but that is not important for our purposes here. 
What is important is that Eqs. (J21))-(I2H1) result in 



ir(°) (g, k]E)^ — for q:>j,k; A; ~ 7. (29) 

g7 

This means that the Eqs. ()19 p - (j^^ are well behaved, in the sense that for A ^ fc, 7 the 

3 

results are independent of A — up to corrections of relative order -^g. This is equivalent to 
the statement that Eq. can be renormalized — i.e. made to render cutoff-independent 
predictions — by the introduction of a momentum-independent force like the one written 
explicitly in Eq. Q. But the set of equations developed here (following Refs. [l3|, l32| ) 
allows us to compute three-body phase shifts without the inclusion of an explicit three-body 
force. 

Furthermore, the energies at which the homogeneous version of Eq. (fTU|) : 

K^^\q,k;-Bs) = f dq' q'' Y^'\q,q'--B,)K^'\q\k--B,) , (30) 

has a non-trivial solution will be the energies at which three-body bound states occur. In 
fact, there are infinitely many such energies, most of which lie outside the domain of validity 
of the EFT, and so should be regarded as spurious solutions of Eq. (IHUj). However, those 
eigenenergies -B3 which are within the radius of convergence of the theory represent predic- 
tions which should be reliable within the error estimate at the order under consideration 
(leading order in this case). The EFT equations developed here therefore explicitly show 
that there will be correlations between the three-body scattering length 03 and bound-state 
energies which are < 1/mB?. 

III. THE SUBTRACTED EQUATIONS AT NNLO 

In this section we will discuss how to include corrections which are higher order in the 
R/a2 expansion using the subtraction formalism. We will show that up to next-to-next-to- 



leading order in this expansion (NNLO) no additional three-body counterterm is needed for 
consistent renormalization. 

Whether a new counterterm is needed or not at a particular order of an EFT can usually 
be determined with the help of naive dimensional analysis (NDA). In perturbative EFTs 
such as chiral perturbation theory this sort of analysis allows us to predict at what order 
terms of a given structure must be included in the Lagrangian in order to renormalize the 
theory. 

However, in EFTs where non-perturbative resummation of a particular set of diagrams 
is mandated by the power counting, NDA may not be a reliable guide to the size of higher- 
dimensional operators. For instance, in the EFT that we are using here, and whose La- 
grangian is given by Eq. Q, NDA leads us to expect that the size of the three-body force 
Dq will be R^/M. This would make that three-body force an effect that only needed to be 
considered at NNLO. However, Bedaque et al. found it necessary to include such a three- 
body force at leading order, so as to render three-body observables cutoff independent jsf. 
In three-body systems with 02 ^ -R renormalization forces the three-nucleon coupling Dq 
to have a larger size than we would expect based on NDA applied at the momentum scale 
l/R. 

The authors of Ref. |^ demanded that with each further order included in the calcula- 
tion the corrections should scale with an additional power of R/a2- Using analytical and 
numerical arguments they concluded that an additional three-body force enters at NNLO. A 
renormalization-group analysis of the perturbations around the limit cycle found in Ref. |3| 
appears to support this view [33[. In contradistinction to the approach of Ref. jj|, here we 
let NDA and the question of renormalizability be the ultimate guides to power counting in 
the three-body problem with short-range interactions. A new structure should be included 
in the effective Lagrangian of the theory if it is necessary for renormalization, i.e. it must be 
included so as to make observables cutoff independent up to the order under consideration, 
or if its presence is implied at that order by NDA with respect to the heavy scales in the 
theory. These criteria have the advantage that they are not subjective. They therefore al- 
low the EFT (j2)) to make testable predictions about data in systems in which there are only 
short-range forces. If, while using the number of short- distance parameters that is indicated 
at a given order by our combined criteria of renormalizability and NDA, the EFT Q is 
unable to describe data with the accuracy expected at that order then we can only conclude 
that the forces at work in that system are not truly short-ranged. 

Analysis in this spirit beyond leading order is most transparent if a perturbative expan- 
sion of the scatterin g am plitude (or, equivalently, the K-matrix) in powers of R/a2 can be 
established. In Ref. 121 Mehen and Hammer showed that in such an expansion no three- 



body force beyond the h term in the leading-order computation of Ref. [3| was needed to 
renormalize the 0{R/a2) piece of the scattering amplitude. 

However, it should also bepossible to imitate the approach of Weinberg in the nuclear 



effective theory with pions J3J, |35| and expand the kernel of the integral equation gov- 
erning 1+2 scattering in powers of Rja^. If we proceed in this way we are certainly not 
doing a calculation that is less accurate than the perturbative one discussed in the previous 
paragraph — provided we can properly renormalize the resulting amplitudes. In fact, here we 
will show that if we calculate the kernel including terms of order i?/a2 and -R^/a| then the 
renormalization already displayed for the leading-order calculation in Section |H] is sufficient 
to facilitate predictions for 1-1-2 scattering. In other words, the same subtraction suffices to 
renormalize the theory at leading, next-to-leading, and next-to-next-to-leading order. 



To develop the integral equation that governs 1+2 scattering at higher orders in the R/a2 
expansion we first note that the full two-body t-matrix with effective-range contributions 
resummed is given by J36|, l37| 



r(E) = -— ; . (31) 

vrm -7 + y/^^^ + I (7^ + mE) ^ ' 

Corrections to Eq. ()31|) involve the shape parameter, and, for natural values of that param- 
eter, are suppressed by [Rja^f" . The r of Eq. (jHTj) has poles which represent bound states 
with energies outside the region of validity of the EFT. It therefore cannot be used in the 
integral equation for cutoffs A > 1/r unless further effort is devoted to the subtraction of 
these unphysical bound states. Therefore, we expand the dimeron propagator (jHTj) up to a 
given order in Rja^. The leading-order result was already displayed in Eq. Q. At higher 
orders we obtain 



4..),^^ _«""(«) 



^""(^) = -wtt' '^^' 



with: 

2 



^^"^i?) = ^E©^[^+v^=^r"^ (33) 



Tim 

■ =0 



although actually at n = 3 and beyond contributions from the shape parameter term, and 
other terms omitted in Eq. (jnU, need to be added in order to get the full result to OilT' jd!^. 
Regardless, once n > 0, the r of Eq. (jH^ no longer has the same residue at i? = — i?2- The 
relationship between the K-matrix and the phase shifts is affected by this change in wave- 
function renormalization. It becomes: 

K^'^Xk, k; E) = -_J-^^ tan5. (34) 

Since K^"\0, 0; — -B2) enters the driving term of our subtracted integral equation for the 
threshold amplitude, not only the kernel, but also the inhomogeneous term, of the integral 
equation, changes as we go to higher order. This, however, is the only change in the driving 
term of that equation, since Z is the same no matter what order in R/a2 we calculate to. 

Now that we have inserted the additional two-body length scale r ^ R into our theory, 
the breakdown scale of the EFT is explicit in our three-body integral equations. In order 
to demonstrate renormalization-group invariance of observables in the theory we must show 
that amplitudes for p <^ 1/r are cutoff independent (up to higher-order terms) for A ^ 1/r. 
Therefore in what follows we analyze the integral equation for the nth-order amplitude in 
the domain p ^ 1/r ^ A. 

Consider the integral equation at threshold at NLO (n = 1) or NNLO {n = 2): 



ir(")(p,0;-52 



2rn, 



p2 _|_ ^2 



4m^ 

+ — — 



with 



d,i|og (4±4±«±4) S<'')(-B,;,)A'«"'ta,0; -B,), 

pq \p + q — pq + "y J 

(35) 
5('^)(-52;g) = 5(")(-52-^), (36) 



where the right-hand side of this equation is given by (joa|) . For p ^ 7, and after a subtrac- 
tion, this equation becomes 



K^^\p,0--B2 



2m 



/yZ 



+ ir(")(0,0;-E2 



+ 



Am 



2 /-A 



dg 



1 , f p +pq + q \ 2 
— log 
pq 



pi _ pq j^ q2 J q2 



S^''\-B,;q)K^^\q,0;-B,). 

(37) 



We must now consider two different momentum domains within the integral on the right- 
hand side of Eq. (jHTj) . In the first we have momenta q <^ 1/r, which are within the domain 
of validity of our EFT, while in the second we have q ^ 1/r, i.e. q is above the EFT 
breakdown scale. To make this transparent we separate the integral in Eq. ()37|): 



K^^\p,0;-B2 



4m2 r^/^ 



ma 



1 , fp^ + qp + q^ 

~ l°g ~2 ^ 

pq \p^ — qp + q^ 



q' 



4m' 



2 /-A 



dg 



1/r 



1 , (p^ + qp + q^\ 2 

~ l°g ~2 ^ ~ ~ 

pz _ qp _^ qz ] nz 



pq 



q^ 



where a is a number of order one: 



3^37 o 
a = 2. 



S^^\-B2;q)K^^\q,0;-B2) 

S^-\-B^-q)K^^\q,Q--B^), 

(38) 

(39) 



Now we take p <C 1/r, giving 
i^(")(p,0;-52) 



ma 4m^ f^^^ , 
T 3 Jo 



1 , fp^ + qp + q^\ 2 

~ l°g ^ ^ ~ ~ 

\p^ — qp + g^/ ^^ 



16m,' 



pg 
,2 



g^ 



For g ^ 1/r we have 



- / dq^^S^'^\-B,;q)K^"\q,0;-B,) 

Ji/r q 



7rm^ \ 4 / 



and we can reduce the second integration to obtain 



K^"\p,0--B2 



m,a 



r 



+ 



4m2 /■!/- 



dg 



1 , fp'^ +pq + q^ 

~ i°s — -^ 

pq \p^ — pq + q'^ 



qz 



VS; 



r 



16 



4 / 3v37r Ji/^ 



^W(-52;g)i^("ng,0;-B2) 

(40) 
(41) 

S^^\-B2-q)K^^\q,Q--B2) 

(42) 



g' 



3— n 



For g ^ 1/r we observe numerically that: 

i^(")(g,0;-i?2) 



r/i 



,/?(") (g) 



7l+n/2 



(43) 
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with P^'^\q) oscillatory and bounded. (See Fig. QJ) Therefore for p <^ 1/r we can write the 
integral equation as 



ir^"^ p, 0; -^2 = ^ + — - / dg — log i^^ — , 2 - - 




3 VStt ii/^ ^'' g^-"/2 



S(")(-i?2;g)i^("n9,0;-i?: 



2) 



_mp'lZ^\ _^ / dg^--^. (44) 



The second integral converges even when the cutoff is taken to infinity. Indeed, for 
A ^ 1/r it is dominated by g ~ 1/r, which allows us to estimate its size to be: 



- „4— n/2 ^ \ r 

1/r y V' 

Therefore, up to corrections of this size, the high-momentum piece of the integral equation 
decouples and 

i^(") {p, 0; -^2) ~ — for 7 < j9 < 1/r, (45) 

IP 

for n = 1 and n = 2, which is the same behavior for p ^ 7 as at leading order. Matching the 
low-momentum behavior of Eq. P3j) to the observed high-momentum behavior ()43|) allows 
us to determine that: 

This finally leads us to conclude that the size of the high-momentum piece {q ^ 1/r) of the 
homogeneous term in the integral equation is ~ —p'^r^, as long as p <^ 1/r. Comparing to 
the behavior of Eq. (jl^ we see that the contribution of high-momentum modes is suppressed 
by three powers of the small parameter pr, and so will only be significant if one wishes to do a 
calculation up to N^LO accuracy. At both next-to-leading and next-to-next-to-leading order 
the threshold amplitude obtained by the subtraction technique will be cutoff independent 
up to terms which are higher order in r/a2 and/or pr. 

In fact, the behavior (j43|) of K^^\q,0; —B2) at large q can be understood via Hilbert- 
Schmidt theory. If we first symmetrize the kernel of our (unsubtracted) integral equation, 
and then make a Hilbert-Schmidt decomposition of the new kernel ^, we find that the solution 
for the threshold amplitude can be written as: 

ir(")(g,g'; -B,) = ^^^\^)^^^^ ^^') _ (47) 

In Eq. (jTTjl {Em : m = 1, . . . , cxd} is the set of eigenfunctions of the integral-equation kernel, 
with each one corresponding to a different eigenvalue Am- These eigenfunctions are only 
normalizable in the A ^ cxo limit if 

/■oo 

/ dqTt\q)S^-\q)Tt\q)<^. (48) 

Jo 



^ Since the procedure for the subtraction is carefully designed to maintain the symmetry of the K-matrix 
the subtraction does not modify the validity of the Hilbert-Schmidt decomposition. 
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FIG. 1: The amplitude i^*^"'^(p,0; —i?2)p"' as a function of the momentum p with r = 0.0737~^, 
as = 1.3627"^ and A = IOOOO7. The sohd, dashed, and dot-dashed Unes denote the LO, NLO, 
NNLO result, respectively. The LO line reproduces the results given in |3l. Il3i|. 

In the case of the homogeneous equation which predicts the eigenenergies at which 1 + 2 
bound states appear the requirement PHj) corresponds to demanding that the bound-state 
wave function is normalizable. 

Assuming that the eigenfunction that is largest at high momentum falls off as T^"^ ~ -^ 

and using the general asymptotic behavior of S at arbitrary order n: S"^"-* ~ q"'~^^, Eq. 
implies: 



n 



n + 1- 23fJ(a„) < -1 ^ 3fJ(a„) > - + 1 



(49) 



Substituting this behavior back into the Hilbert-Schmidt decomposition of K^"'\ Eq. (07j), 
we see that K^"'\q, q'; —B2) must fall off at least as fast as ^J^+i for l\ 1 ^ Q- Numerically 
we find that the equality in Eq. ()49p is realized at LO, NLO, and NNLO. But as we have 
seen previously, such a fall-off with off-shell momentum q is already enough to guarantee the 
decoupling of the high-momentum {q ^ 1/r) modes in the integral equation that governs 
the threshold scattering amplitude. 

We now turn our attention to the amplitude when E > —B2. The NLO and NNLO 
versions of Eqs. (|T^ - P^ are found by making by the replacements: 



K(o)^ 



q,q 



T 



(0) 



-B2) 

{E;q) 



K^^\q.q"-. 



-Bo 



T 



(50) 

(51) 



The question is whether these equations yield cutoff-independent predictions for the ampli- 
tude K{p, q; E) when p,q,k <t^ 1/r. 

Since the fall-off of 6Z at large momenta is, at worst, that given in Eq. (j28j) we see that 
as g — *> cx), 6Z{q',q; E) falls off at least as fast as Tm (q) — at least for n = 1 and n = 2, 
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which are the only cases we consider here. In consequence, as long as the integral: 

jyq'q"Tt'\qV-\E;q')Tfiq') (52) 

is not cutoff-dependent then Eqs. |T9|) - p3|) will yield cutoff-independent results. But we 
already know from Eq. (j48p that this integral is convergent in the A — >■ cxd limit, and so we 
are guaranteed that, at least at NLO and NNLO, the formalism of the previous section can 
be adapted so that, at any energy such that k <^ 1/r, it yields cutoff-independent results 
for the low-momentum amplitude. 

We note, in passing, that the integral (j3^ would not be convergent if Tm (q) ~ l/q as 
q -^ oo, i.e. the NLO and NNLO amplitudes had the same fall-off at g ^ 1/r as the LO 
amplitude. If the LO high-momentum behavior were also present in the higher-order cases 
then the integral of Eq. ()52|) would behave as r^/c^lnA for large A. The low-momentum 
amplitude would then not be renormalization-group invariant, and a counterterm, like the 
one introduced in Ref. jj], would be necessary to remove this divergence. However, our 
numerical results (see Section I1V|) indicate no A-dependence (to six significant figures) in 
the NNLO phase shifts, because the fall-off of K'-^^q, 0; —B2) with q is really l/g^+"/^ and 
not 1/g. 

Therefore we have shown that with the changes in S written in Eq. ()33|) . and the cor- 
responding changes in the inhomogeneous term, we can use Eqs. fll5|) and (fT ^ -(|^ ^ to 
compute three-body scattering and bound-state observables at NLO and NNLO. In these 
equations we can take the cutoff to be arbitrarily large. Therefore we have proven — using 
a combination of numerical and analytic arguments — that no new counterterm is needed 
to renormalize the three-body system with short-range interactions at NNLO. This conclu- 
sion is driven by the behavior of the zero-energy scattering amplitude for off-shell momenta 
q ^ 1/r, and that behavior is tied to the requirement that the bound-state wave functions 
of the 1 + 2 system be normalizable. 

IV. THE THREE-BOSON SYSTEM AT NNLO 

As mentioned above, few-body systems of ^He atoms seem to be an ideal testing ground 
for the formalism laid out here. Helium-4 clusters are bound by the van der Waals force, 
a potential with a long-range tail of the form —Cq/t^. The S-wave phase shifts for this 
potential obey a modified effective-range expansion: k cot 6 can be expanded in powers of 
fcf, but non-polynomial terms appear at order k^. As pointed out by Braaten and Hammer |9|, 
38|, this means that an EFT with contact interactions should be able to mimic the underlying 



van der Waals interaction up to this order, even though the non-analytic effects at O(fc^) 
cannot be reproduced by a short-range potential. The size of the 0{k'^) term in the modified 
effective-range expansion is / ~ {mCQ/h'^Y^^, and so, when applied to these systems, the 
EFT with contact interactions is an expansion in l/a2 and kl instead of R/a2 and kR. In 
the following we will present numerical results for the "^He three-body system up to NNLO 
in this expansion. 

The atom-atom scattering length 02 = (l04l^g) A and the dimer binding energy B2 = 
1.30962 mK were derived from measurements of the dimer bond length, and using the zero- 
range approximation [39| this gives an effective range of the order of 10 A. Therefore the 
EFT expansion parameter is l/a2 ~ 0.1 and the EFT expansion should converge rather 
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Input 




b^'\b,] 


b'^\b,] 


«3[7-^] 


-°3 


LO 
NLO 

NNLO 


1.738 
1.738 
1.738 


99.27 

84.87 
89.52 


1.179 
1.199 
1.203 


03 


LO 

NLO 

NNLO 


1.723 
1.735 
1.737 


97.12 
84.48 
89.38 


1.205 
1.205 
1.205 


TTY[18,J^ 




1.738 


96.33 


1.205 



TABLE L EFT predictions for the ^He trimer binding energies and atom-dimer scattering length 
as up to N^LO. Energies and lengths are given in units of two-body binding energy. The first three 
rows show LO, NLO, and NNLO results when B^ is used as input. The next three rows give the 
results when 03 is chosen as input. In the last row we display the results of the calculations of 
Ref. BH- 



Input 




B!^hB2] 


B^m 


ash-'] 


B^l' 


LO 
NLO 

NNLO 


1.741 
1.741 
1.741 


99.78 
85.31 
90.08 


1.173 
1.193 
1.196 


as 


LO 
NLO 

NNLO 


1.643 
1.655 
1.656 


85.02 
74.79 

78.82 


1.362 
1.362 
1.362 


TTY[20] 




1.741 


96.06 


1.362 



TABLE II: EFT predictions for the He trimer binding energies and atom-dimer scattering length 
03 up to N^LO. Energies and lengths are given in units of two-body binding energy. The first three 
rows show LO, NLO, and NNLO results when B^ is used as input. The next three rows give the 
results when 03 is chosen as input. In the last row we display the results of the calculations of 
Ref. ET 



quickly. Since we compute observables up to NNLO we expect our results at that order have 
a remaining error ~ 0.1%. 

To compute observables in the three-body system we need one three-body datum. How- 
ever, the experimental information on clusters of ^He bound states with more than two 
atoms is rather limited. While the '^He trimer and larger clusters have been observed, no 
quantitative information about their binding energies is available llSl Il6l. On the other 
hand, various "realistic" potentials have been developed (see, e.g. J22,|23|J) and have been 
employed for few-body calculations |l7|, |l8|, |l9|, |20|, |2l|. In need of a three-body observ- 
able to fix the value of our three-body parameter we will employ these calculations as our 
three-body input. 

This strategy follows that of Braaten and Hammer '3^ , who analyzed universal properties 
of three-body systems of ^He atoms using the EFT with contact interactions alone. They 
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computed various observables, including the trimer energies and the atom-dimer scattering 
lengths, using the results obtained with different atom-atom potentials as input parameters. 
They achieved the desired accuracy of 10% for a leading-order EFT calculation and estimated 
the NLO corrections for low-energy observables by assuming they would shift the leading- 
order result by a factor proportional to l/a2. 

The "TTY potential" 22| is one of the potential models used to predict the binding 



energies of clusters of ^He atoms. As with all such models, it relies on theoretical assumptions 
about the underlying atom-atom interactions. It reproduces the experimental value of the 
dimer binding energy given above. The effective range has been calculated for various 
potentials by Janzen and Aziz, and for the TTY potential they obtain r = 7.329 A |40|. 
When used in a three-body calculation the TTY potential produces exactly two three-body 
bound states. There have been several calculations of these trimer binding energies, B^ 
and B\ , employing the TTY potential, and all of them agree on those energies to within 
0.5%. 

The situation for the atom-dimer scattering length is somewhat less clear. Two recent 
Faddeev calculations which both employ the TTY potential a gree on the trimer binding 
energies, but have values for 03 which disagree at the 10% level 19|,|20|. This is more than 
the combined errors stated in Ref. (5 A?s 4%) and Ref. \l^ fO.l Kk. 0.1%). In what 
follows we will examine both Faddeev calculations in light of the EFT treatment developed 
in the previous sections. Although we focus largely on the results of these calculations when 
the TTY potential is their two-body input, the pattern of results that we find is the same if 
we compare to calculations that employ other atom-atom potentials instead. 

Our results for the TTY potential using input from the computation of Refs. |l8|, |l9| 
are shown in Table Ql Note that for better comparison we have computed all bound state 
observables in units of the two-body binding energy B2. At low orders the results change 
somewhat depending on whether we use the atom-dimer scattering length or the energy of 
the higher-lying trimer Bj^ as the three-body input datum, and so in the table we present 
the numbers we obtained via both strategies. In this case the dependence on the choice of 
input becomes systematically less as the calculation is carried to higher orders — one sign 
that the EFT is converging in a self-consistent way. Our computations are accurate to at 
least the number of digits shown in the table. The results are also cutoff independent for 
sufficiently large cutoffs, which supports the claim of the previous section that no additional 
three-body counterterm is needed for the renormalization of the three-boson problem at 
NNLO. 

At NNLO we achieve very good agreement for the energy of the excited three-body bound 
state or the atom-dimer scattering length if we use (respectively) 03 or B)^ as calculated 
in jl9|]. If we use S3 as the input parameter we obtain 03 = 1.2037^^, to be compared 
to 1.2057"^ obtained by Roudnev. The remaining discrepancy — which is significant given 
the accuracy of our calculation and of that of Ref. |l9| — is about 0.1%, exactly what we 
expected from a calculation at this order. The situation for the deep-lying trimer is less 
pleasing. We obtain 89.38i?2 for the three-body ground state, which has to be compared to 
Roudnev's full Faddeev TTY-potential result of 96.33i?2- The disagreement between these 
results for the lowest-lying three-body bound state appears to be due to slow convergence of 
the EFT for that observable. This binding energy changes by roughly 15% from LO to NLO 
and 5% from NLO to NNLO. (The fact that the NLO correction is large, and moves the 
EFT result away from the calculation of Ref. ;20] has also been found in an EFT calculation 

that uses the strictly perturbative l/a2 expansion |41.].) Since S3 is large on the scale of B2 
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FIG. 2: The Phillips line for leading (dot-dashed line), next-to- leading (dashed line) and next- 
to-next-to-leading (solid line) order. The cross and diamond correspond to the TTY and HFD-B 
prediction obtained in llgl. |l9|, respectively. The square and triangle denote the corresponding 
predictions obtained in |20|. 



this relatively slow convergence is not surprising, although the pattern of NLO and NNLO 
corrections suggests that an N^LO calculation of B^ should have at worst a 2% error. 

Since Roudnev & Yakovlev and Motovilov et al. essentially agree on the value of -B3 , 
but disagree on the number for 03, the fact that we reproduce the 03 result of Ref. [19| means 
that we cannot also reproduce the numbers of Motovilov et al.. This can be seen explicitly 
in Table |n] If we use, for example, their result for S3 as our three-body input we obtain 
1.1967"^ for the atom-dimer scattering length, as compared to their 1.3627~^. In contrast 
to the result for B^ the EFT result for 03 shown in Table ITTl is converged at NNLO, but it 



does not reproduce the number quoted in Ref. |20|. We see this as a strong indication that 
the atom-dimer scattering lengths obtained by Motovilov et al. have a larger error than the 
5 A quoted by the authors. 

A further way to illustrate this point is to consider the Phillips line. The Phillips line 
is a universal feature in three-body systems with large two-body scattering lengths: it is 
a nearly linear correlation between the 1-1-2 scattering length and the three-body binding 
energy. In Fig. |21 we display the change of the Phillips line with the inclusion of higher- 
order corrections. The Phillips line is converging to a definite result for a particular value 
of r/a2. In Fig. |21 we also display the data points for the TTY potential and the HFD-B 
potential as given in J20| and |l8|, |l9| ^. It is clear that the values obtained by Roudnev 



3 The effective range of the HFD-B potential is 7.28 A ^ and therefore differs by 1% from the TTY 
effective range. However, this difference amounts to a 0.1% shift in observables which is negligible on the 
scale of this plot. 
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FIG. 3: The LO (dot-dashed hne), NLO (dashed hne), and NNLO (sohd hne) S-wave scattering 
phase shifts Sq using the TTY results for B2 and r and the TTY prediction of Ref. [20| for the 
atom-dimer scattering length. The filled circles show the results for 6o{E) given in |2C 



and Yakovlev lie very close to the line, while the values obtained by Motovilov et al. seem 
to be systematically shifted away from it. This appears to support Roudnev's argument in 
Ref. jlSl] that the 03 values of Motovilov et al. were not calculated on a large enough spatial 
grid. Very accurate phase-shift calculations at energies close to the atom-dimer threshold 
are required for a precise determination of 03. 

We also computed the phaseshift 5q with 03 = 1.867"^ as the three-body input for our 
subtraction. In Fig. El we display our results for LO, NLO, and NNLO in each case and 
compare it to the TTY results given in Refs. [20|. In this figure we display Motovilov et aUs 
phase shifts because the criticism leveled by Roudnev at their calculation should become 
less of an issue as the energy increases away from the two-body threshold. 

There are a few interesting points to note regarding Fig. |H1 First, by NNLO the EFT 
seems to have converged to a definite result for 6o{E) in this range — once 03 is fixed. Second, 
the NNLO EFT results agree well with the results of the Faddeev calculations using the 
phenomenological TTY potential. The use of the atom-dimer scattering length obtained 
in J20| as the three-body input accounts for whatever difficulties they may have had in 
computing phase shifts close to two-body threshold. Third, the NLO and NNLO corrections 
definitely make the energy-dependence of the EFT results significantly closer to the energy- 
dependence of the TTY-potential phase shifts. Since the slope of the phase shift at zero 
energy is fixed to be 03, and so is independent of the order of the calculation, it is not 
surprising that the effect of the higher-order corrections is most noticeable at the higher 
energies shown. 
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V. SUMMARY 

In this work we have extended the subtraction formahsm developed previously [13|, |3^ 
to higher orders in the R/a2 expansion. Using analytical and numerical arguments we have 
shown that the three-body system at NNLO can be renormalized without the need for 
a second three-body datum. Therefore, one does not need to include an energy- dependent 
three-body force at this order if one accepts renormalizability and naive dimensional analysis 
as the guidelines for the establishment of a consistent power counting. This is an interesting 
result as it provides another example of the potential predictive power of the EFT with 
contact interactions alone. 

Few-body systems of spinless bosons like ^He atoms are the simplest testing ground for 
this EFT. We have advanced the leading-order EFT calculation by Braaten and Hammer 
by two orders — to NNLO — and found interesting results. 

Without a three-body input at leading order the EFT with contact interactions alone is 
not able to deliver any information about the three-body system. Thus, we used a three- 
body observable calculated with a specific phenomenological inter-atomic potential ("the 
TTY potential") to fix our three-body input. Using two two-body inputs and one three- 
body input we achieved very good agreement with the Faddeev calculations of Roudnev 
and Yakovlev [l8|, |l9| for the TTY potential. The disagreement between our calculation 
and the results obtained by Motovilov el. al is of the order of 10%, which is two orders 
of magnitude larger than the expected accuracy of an NNLO calculation. Our calculation 
therefore supports the claim of Ref. [19| that there were problems with the numbers quoted 
in Ref. J20|, and that the error given there should be enlarged. We believe that this is an 
example of EFT's ability to check the consistency of different phenomenological calculations. 
The numbers of Ref. ^] for a^ are simply not consistent with the model-independent pattern 
of r/a2 corrections to the Phillips line that we have calculated using the short-range EFT. 

Furthermore, it is worth noting that although various inter-atomic potentials give similar 
results for low-energy observables in the two-body system, they give different results for 
three-body observables — albeit in correlated ways that are constrained by relations like 
that shown in Fig. |21 But, to decide which two-body potential gives the most appropriate 
description of the three-body system, an experimental measurement of at least one three- 
body observable is needed. The measurement of one such number (e.g. B^ ') would therefore 
be a valuable test of the current realistic potentials and would further our understanding of 
the atom-atom interaction. The results of this paper show that once one knows, for instance, 
i?3 , the EFT can be used to generate precise, model-independent, predictions for remaining 
low-energy three-body observables. 

A further application of the formalism presented in this work is the three-nucleon sector. 
The nucleon-nucleon force is of finite range and so there the effective-range expansion is 
analytic to all orders in fc^. We have seen here that the short-range EFT converges well 
for systems of three Helium-4 atoms, and it has been shown that this EFT also converges 
as expected in the two-nucleon sector, where some observables have been computed up to 
four orders in the R/a2 expansion [2J, |25|. The extension of the equations developed here 



to the the three-nucleon problem requires only the addition of spin and isospin degrees of 
freedom [13|, |4^. This will open the way to EFT computations of low-energy three-body 
observables (including those involving electroweak probes) which have errors of < 5%. 
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